Experimental evidence of the ferroelectric phase transition near the A— point in liquid 

water. 
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We studied dielectric properties of nano-sized liquid water samples confined in polymerized sili- 
cates MCM-4f characterized by the porous sizes ~ 3 — Wnm. We report the direct measurements 
of the dielectric constant by the dielectric spectroscopy method at frequencies 25Hz — 1MHz and 
demonstrate clear signatures of the second-order phase transition of ferroelectric nature at tem- 
peratures next to the A— point in the bulk supercooled water. The presented results support the 
previously developed polar liquid phenomenology and hence establish its applicability to model 
actual phenomena in liquid water. 



Large dipole moments of individual molecules and the 
high density of water are reasons behind the very rich 
phase diagram and quite a few "anomalous" properties 
P~H7]. Understanding the dielectric response of the liquid 
plays a crucial role in modeling molecular interactions in 
computational physical chemistry, biophysics, and drug 
design applications [5J[S]. The static dielectric constant e 
of water at room temperatures is very large, increases as 
the temperature decreases, and even diverges in super- 
cooled water being extrapolated to the unreachable tem- 
perature Tc ~ 228° K of the so called A-transition |10j . 
Most other thermodynamic quantities, such as isother- 
mal compressibility, density, diffusion coefficient, and vis- 
cosity, are also singular [TUl - fT3] . Following the earlier 
idea of [T5], the authors of [TS1 [TB] indicated that the 
phase transition could have ferroelectric features. The 
ferroelectric hypothesis was also supported by a number 
of molecular dynamics (MD) simulations flTT - fliT . For 
example, a ferroelectric liquid phase was observed in a 
model of the so called "soft spheres" with static dipole 
moments |17[ I20H24] . The conclusion seems to be model 
independently confirmed in the MD of hard spheres with 
point dipoles |20l 125] , soft spheres with extended dipoles 
[26] . and in a model of fluids made of two-state particles 
with a nonzero dipole in the excited state |27| . 

Unfortunately, the relation of the ferroelectric phase 
transition to the A— point in actual water or even the 
very existence of the paraelectric phase may be difficult to 
confirm in MD simulations with finite number of atoms. 
There are quite a few reasons to it. First of all due to 
the long range nature of the dipole-dipole interactions 
between the molecules, the simulated liquid tend to form 
tightly correlated domains and the calculated properties 
depend strongly on the boundary conditions in any rea- 
sonably sized system [53]. In the same time because of 
the low temperatures and the proximity of various phase 
transitions the relaxation processes take a very long time, 
which transforms into a necessity to simulate large and 
strongly interacting molecular systems in realistic force- 
fields for a very long time. 

The same problems naturally hinder theoretical un- 
derstanding of the physics behind the phase transition. 



There could be no purely electrostatic MD model of 
the phase transition, since classical systems with electro- 
static interactions only are inherently unstable. Chem- 
ical forces, such as hydrogen bonds, are known to play 
a very important role in water molecule ordering at all 
temperatures in general, but particularly when close to 
the phase transition [281 [25] , The "minimal" continu- 
ous model capable of predicting finer effects depending 
both on the hydrogen bonding properties and the electro- 
static interactions of the water molecules was proposed in 
|30H32| . One of its predictions is the ferroelectric phase 
transition in liquid water within the temperature inter- 
val Tc = —37 -= — 47°C, which is remarkably close to 
the A— point. In [TU] the singularity of dielectric con- 
stant e as extrapolated to the unreachable temperature 
Tq = — 45°C was indeed reported. On the other hand 
due to its weak character and unattainability of Tc the 
authors of [TU] did not attribute it to the a ferroelectric 
phase transition. 

Fortunately the freezing temperature can be essentially 
lowered in water confined in nanopores (see EJ 1331 134] 
and references therein). The porous systems can be 
studied by vapor pressure measurements [35J, dilatom- 
etry [35] |37], calorimetry [36-42 , nuclear magnetic res- 
onance [35] G3S HUSH], neutron diffraction gU H5H5J], 
small-angle neutron scattering |48II53|. and X-ray diffrac- 
tion [5J . The measured transition temperature depression 
can be approximated as AT ss —K/(R — t), where R 
is the characteristic radius of the pore, K = 52K-nm, 
and t — OAnm is the (empirical) hydration layer thick- 
ness [33]. Already for the pore radii below R ~ 2nm 
the decrease in freezing temperature can be as large as 
— and hence the predicted ordering phase transition 
at T = Tc can be observed and studied experimentally. 
Moreover, in hindsight we believe that the first signatures 
of such a transition were already observed in [5] . In this 
work we present the results of the direct measurements 
of the low-frequency liquid water dielectric constant, e, 
in nanopores by dielectric spectroscopy method. The ob- 
served singularity of e near the A— point temperature 
in the supercooled bulk water is a clear signatures of 
the second-order transition of ferroelectric nature, both 
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in quantitative and qualitative agreement with the polar 
liquid phenomenology. 

The polar liquid phenomenology [21] extends the con- 
tinuous models [211 I54TI6T] and was originally developed 
to describe electrostatic energies of biomolecules in aque- 
ous solutions for drug discovery applications [62J. The 
model naturally describes the ordering phase transition 
in water |30l I3"2"] . Within the suggested model the po- 
lar liquid is characterized by the vector-field s(r) = 
(d (r)) /do, where d (r) is the vector of the static dipole 
moment of molecule residing at point r, do is its absolute 
value, < s(r) < 1. Total dipole moment of molecule 
equals d f (r) = d (r) + d e (r), where d e (r) is the dipole 
moment induced in the electronic shell. The Helmholtz 
free energy of polar liquid is described by the functional 
of independent variables s (r) and d e (r): 



F(s(r),d e {v)) = P 2 J dV 
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Here C is the phenomenological parameter responsible 
for the iJ-bond network rigidity, Po — nodo, no is the 
particle density of the liquid, £ m is the part of the 
dielectric constant of the liquid which is unrelated to 
the molecular rotations and accounts properly for the 
polarization of the internal (e.g. electronic) degrees 
of freedom. It is the dielectric constant at frequen- 
cies larger than that correspondent to reorientation of 
molecules. Next, Ep = — V</?p is the polarization elec- 
tric field produced by the polarization charges with a 
density pp — —VP (sometimes it is called as the de- 
polarizing field), E e (r) = — V<^ e is the external electric 
field induced by external charges with the density p e (r) , 
P (r) = Po s ( r ) + "-od e (r) is the polarization vector of 
liquid at point r, E (r) = Ep (r) + E e (r) = — Vif is the 
total electric field at point r, <p = tpp + Lp e . Electric po- 
tentials ipp and ip e should be found from correspondent 
Poisson equations Atpp — —Anpp, Aip e = — 47rp e . 

The phenomenological dimensionless function V(s 2 ) 
serves as the polar liquid equation of state and can not be 
established in general form. Its specific form can be found 
by comparing the model results obtained with the help of 
Eq.Q with the results of MD simulations of a specific liq- 
uid. On the contrary, in the small s 2 < 1 limit the func- 
tion takes a nearly universal form: V(s 2 ) ~ As 2 /2 + Bs A , 
A = 47rT/(3eoo) has a universal (i.e. the same for any po- 
lar liquid) form, B = 37r /(5eoo) is a liquid-specific con- 
stant, B ~ 1, and r = (T — T c )/T c , where, at last, 
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is the critical temperature within the model. The phys- 
ical picture behind the change of the coefficient A sign 
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Figure 1: Temperature dependence of low-frequency real part 
of dielectric permittivity of liquid water (solid curve) and ice 
(dashed) in MCM-41 pores with diameters 3.5nm. 



can be recovered from the following argument. For uni- 
formly polarized liquid, s(r) = const, the free energy Q 
takes Landau-like form 163 : 



F = V^( 
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At temperatures T > Tq the equilibrium state corre- 
spond to the disordered paraelectric phase with (s) = 0, 
whereas at lower temperatures T < Tc the polar liquid 
undergoes the second order phase transition and trans- 
forms to the long-range ordered ferroelectric state. 

Second order ferroelectric phase transition should man- 
ifest itself as a singularity in the liquid dielectric constant 
e. To analyze the dielectric response we apply a uniform 
electric field E and calculated the polarization of the liq- 
uid by the minimization of free energy : 



2V (s 2 ) s = — E, d e 
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Consequently, in the weak electric field limit E -c Pq at 
|t| <C 1 the static dielectric constant of the liquid is given 
by: 

l + ^j,T>T c ;e = e OQ (l+^j, T <T c . 

Therefore the measurements of the temperature depen- 
dence in e(T) should exhibit a A-point feature and diverge 
at T = Tq- The exact value of the critical temperature 
for water can be obtained using either of a few published 
measurements of the asymptotic values: e x m 4.9 from 
|64| . eoo = 5.1 from jHS], or = 5.5 from Respec- 
tively, Eq. (|) gives T c = 236K (-37°C), T c = 226K 
(-47°C), and T c = 210K (-63°C). All the numbers are 
remarkably close to T c w 228K [TMT2] or T c « 23 IK 
measured in supercooled bulk water |13| . 

To resolve the issue we investigated the dielectric re- 
sponse of the water samples confined in polymerized sil- 
icate MCM-41 characterized by a typical pore diameters 



3 



of D ~ 3 4- lOnm |67l 155] . We used the dielectric spec- 
troscopy method at frequencies range 25Hz-lMHz. The 
results of the measurements are summarized on Fig[I]and 
show a very distinguished A- feature at T^ xp « —35^x7 in 
full accordance the second-order phase transition picture 
presented above. At first glance at large temperatures 
the data seem to suggest a rise of e with temperature 
in contradiction to both earlier experiments [10 and the 
theoretical prediction Q. In fact the measurement er- 
rors in the dielectric spectroscopy experiment are quite 
large and the contradiction disappears entirely when the 
close vicinity to the transition point is considered. Next 
to the observed transition temperature the singularity of 
e is much stronger stronger than the previously reported 
dependence e oc |r| _a characterized by the critical index 
a w 0.13 [TU]. The mean field theoretical prediction Q 
a = 1 results in a much better agreement with the mea- 
sured values. Close to the phase transition one can con- 
sider using a more more refined approach. Indeed, when 
\t\ <C 1, the fluctuations in the model (JTJ are "force"-less, 
V • s = 0, pp = 0, Ep = |31) . the dipole-dipole interac- 
tions vanish and the scale invariant calculation from [55] 
gives a ~ 1 + (4 — d) /6, where d is the number spatial 
dimensions. Actually d = 3 and a = 7/6 ~ 1.2. which 
gives is stronger singularity than in the mean field the- 
ory. Though it is hard to tell from a few points we have, 
the experimental divergence of e on Fig. [T] appears even 
stronger than each of the theoretical predictions. This 
should not be surprising since our model is indeed over- 
simplified. The calculations can be further improved by 
including at least the two important forms of the liquid 
water states (such as the hexagonal and the cubic water 
structures as explained in e.g. |70[ ITT]). If we associate 
the ferroelectric phase transition with a single (say, cu- 
bic |16| ) component only, then the critical behavior of 
the dielectric constant may be changed quite dramati- 
cally by the sharp temperature dependence of the cubic 
water fraction |72] , 

The hydration thickness Rd <C t < Lt, where Rd ~ 



0.2nm, and Lt ~ 0.7nm are characteristic phcnomcno- 
logical scales of the liquid characterizing roughly the wa- 
ter molecules cluster and the maximum correlated do- 
main sizes respectively [30 32 . Thus the pores in MCM- 
41 are very large compared to the hydration layer thick- 
ness, D 3> t,Lx, and therefore the water in nanopores 
still has all the properties of the bulk water. The equilib- 
rium freezing temperature of water in MCM-41 systems 
with D = 3.5nm is -38°C < Tg* [33] and therefore the 
water samples in our experiments are still liquid at tem- 
peratures T ~ Tq p near the A-point. The liquid and the 
solid water states were definitely distinguished using the 
hysteresis effect EE1 E3] : the solid curve on Fig. l] cor- 
responds to the supercooled liquid state and the dashed 
line describes the dielectric properties of the overheated 
ice state. 

In this way water freezing temperature suppression and 
the hysteresis phenomena typical for water in nanopores 
let us actually investigate the properties of the super- 
cooled bulk liquid water in the temperature range —704 
+23°C. The dielectric spectroscopy measurements reveal 
the strong singularity in dielectric constant at tempera- 
tures corresponding to the position of A-point of the bulk 
water. Therefore the experimental data point support 
the earlier theoretical predictions |30| 152] of existence of 
the second order phase transition of ferroelectric nature. 
Eq.Q reconcile earlier theoretical predictions of (23] with 
the experimentally observed values for the phase transi- 
tion temperature by introducing electronic polarizations 
too. The very fact that the dielectric constant singular- 
ity was actually observed suggests that experiments with 
water in nanoporous materials can be used for in-depth 
studies of complicated phase diagram of water. 
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